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Abstract 

We investigate quantum tunnelling methods for calculating black hole 
temperature, specifically the null geodesic method of Parikh and Wilczek 
and the Hamilton-Jacobi Ansatz method of Angheben et al. We con- 
sider application of these methods to a broad class of spacetimes with 
event horizons, inlcuding Rindler and non-static spacetimes such as Kerr- 
Newman and Taub-NUT. We obtain a general form for the temperature 
of Taub-NUT- Ads black holes that is commensurate with other methods. 
We examine the limitations of these methods for extremal black holes, 
taking the extremal Reissner-Nordstrom spacetime as a case in point. 

1 Introduction 

There are several methods for deriving Hawking radiation [I]-jl5) and for calcu- 
lating its temperature. The original method considered the creation of a black 
hole in the context of a collapse geometry, calculating the Bogoliubov transfor- 
mations between the initial and final states of incoming and outgoing radiation. 
The more popular method of analytic continuation to a Euclidean section (the 
Wick Rotation method) emerged soon after. Relying on the methods of finite- 
temperature quantum field theory, an analytic continuation t ~* it oi the black 
hole metric is performed and the periodicity of t (denoted by f3) is chosen in 
order to remove a conical singularity that would otherwise be present at fixed 
points of the U{1) isometry generated by d/dr (the event horizon in the original 
Lorentzian section). The black hole is then considered to be in equilibrium with 
a scalar field that has inverse temperature (3 at infinity. 

Recently a semi-classical method of modeling Hawking radiation as a tun- 
neling effect was proposed inj-Cni- This method involves calculating the imagi- 
nary part of the action for the (classically forbidden) process of s-wave emission 
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across the horizon (first considered by Kraus and Wilczek 0-13), which in turn 
is related to the Boltzmann factor for emission at the Hawking temperature. 
Using the WKB approximation the tunneUng probability for the classically for- 
bidden trajectory of the s-wave coming from inside to outside the horizon is 
given by: 

r cx exp(-2Im/) (1) 

where I is the classical action of the trajectory to leading order in H (here 
set equal to unity) Expanding the action in terms of the particle energy, 
the Hawking temperature is recovered at linear order. In other words for 
2I = /3E + 0{E^) this gives 

r ~ exp(-2/) ~ exp{-PE) (2) 

which is the regular Boltzmann factor for a particle of energy E where P is the 
inverse temperature of the horizon. The higher order terms are a self-interaction 
effect resulting from energy conservation [01,11]; however, for calculating the 
temperature, expansion to linear order is all that is required. Two different 
methods have been employed to calculate the imaginary part of the action - one 
used by Parikh and Wilczek |5] and the other by Angheben, Nadalini, Vanzo, 
and Zerbini (14) (which is an extension from the method used by Srinivasan and 
Padmanabhan |13p. 

The former method considers a null s-wave emitted from the black hole. 
Based on previous work analyzing the full action in detail the only part 

of the action that contributes an imaginary term is J^"^'' prdr, where pr-is the 
momentum of the emitted null s-wave. Then by using Hamilton's equation and 
knowledge of the null geodesies it is possible to calculate the imaginary part of 
the action. We will refer to this approach as the null geodesic method. 

The latter method involves consideration of an emitted scalar particle, ig- 
noring its self-gravitation, and assumes that its action satisfies the relativistic 
Hamilton-Jacobi equation. From the symmetries of the metric one picks an ap- 
propriate ansatz for the form of the action. We will refer to this method as the 
Hamilton-Jacobi ansatz. 

In this paper we examine these two methods in the context of a broader 
class of spacetimes than has previously been studied. One of our prime motiva- 
tions is to understand the applicability of the method to stationary black hole 
spacetimes such as the Kerr-Newman and Taub-NUT spacetimes. The Taub- 
NUT metric is a generalization of the Schwarzschild metric and has played an 
important role in the conceptual development of general relativity and in the 
construction of brane solutions in string theory and M-theory. jJS] The NUT 
charge plays the role of a magnetic mass, inducing a topology in the Euclidean 
section at infinity that is a Hopf fibration of a circle over a 2-sphere. "A counter 
example to almost anything" , Taub-NUT spaces have been of particular 
interest in recent years because of the role they play in furthering our under- 
standing of the AdS-CFT correspondence |^,|^,|2nj. Along these lines, the 
thermodynamics of various Taub-NUT solutions has been a subject of intense 
study in recent years. Their entropy is not proportional to the area of the event 
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horizon and their free energy can sometimes be negative JSIjUOI) |22|. 
Solutions of Einstein equations with a negative cosmological constant A and 
a nonvanishing NUT charge have a boundary metric that has closed timelike 
curves. The behavior of quantum field theory is significantly different in such 
spaces, and it is of interest to understand how ADS/CFT works in these sorts 
of cases [33] . 

All such thermodynamic calculations have thus far been carried out in the 
Euclidean section, using Wick rotation methods. For most Taub-NUT spaces 
the Lorentzian section has closed timelike curves. As a consequence, determina- 
tion of the temperature via the original method of Hawking - while mathemat- 
ically clear - is somewhat problematic in terms of its physical interpretation. It 
is straightforward enough to analytically continue the time coordinate and var- 
ious metric parameters to render the metric Euclidean. Regularity arguments 
then yield a periodicity for the time coordinate that can then be interpreted as 
a temperature. However the Lorentzian analogue of this procedure is less than 
clear, though it has been established that a relationship between distinct ana- 
lytic continuation methods exists |21]. An independent method of computing 
the temperature associated with event horizons in NUT-charged spacetimes is 
certainly desirable. 

Our goal in this paper is to address this question, and to more generally inves- 
tigate the tunnelling approach outside of the spherically symmetric ansatz. To 
this end, we compare the null geodesic method and the Hamilton- Jacobi ansatz 
for obtaining the imaginary part of the action. We then apply these methods to 
a variety of spacetimes, and derive a general formula for the temperature from 
this method. We then consider specific cases of interest, beginning with Rindler 
space and moving on to charged and rotating black hole spacetimes. Turning to 
Taub-NUT spaces, we obtain a general expression for the temperature for a sub- 
class of Taub-NUT spacetimes without closed timelike curves (CTCs) that we 
can compare to those obtained via Wick rotation methods. We find agreement 
in all relevant cases. 

Our paper is structured as follows. The next section will outline the two 
methods, starting with a discussion and generalization of the null geodesic 
method and followed by a discussion of the Hamilton- Jacobi ansatz. We demon- 
strate that knowledge of the total mass or energy is not essential by showing 
the direct application of these methods to Rindler spacetime. We then apply 
these methods to stationary space-times, considering in turn the Kerr-Newman 
class of metrics and then Taub-NUT spacetimes. In each case we obtain re- 
sults commensurate with other methods, concentrating in the latter case on the 
subclass of Taub-NUT- AdS spacetimes that do not have closed timelike curves 
|23| . We finish with a preliminary discussion of issues that occur when applying 
the method to extremal black holes, concentrating on the specific case of the 
extremal Reissner-Nordstrom spacetime. 
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2 Calculating the Imaginary Part of the Action 
For an Outgoing S-Wave 

2.1 Null Geodesic Method 

We begin by reviewing the null geodesic method The general static spher- 
ically metric can be written in the form 

ds^ - -f{r)dt^ + ^ + r^dn^ (3) 
9{r) 

which covers a broad range of black hole metrics. We want to write it in Painleve 
form ^3 so that there is no singularity at the horizon. This is easily accom- 
plished via the transformation 

' — //Ml* 

yielding 

ds^ = -f{r)di^ + 2y/J{r)J-^ - Idrdt + dr^ + r^dU^ (5) 



This coordinate system has a number of interesting features. At any fixed 
time the spatial geometry is flat. At any fixed radius the boundary geometry 
is the same as that of the metric (jSJ . 

The radial null geodesies for this metric correspond to 



±l-^T^W)] (6) 



V 9{r) 

where the plus/minus signs correspond to outgoing/ingoing null geodesies. 

The basic idea behind this approach is to regard Hawking radiation as a 
quantum tunnelling process. However unlike other tunnelling processes in which 
two separated classical turning points are joined by a trajectory in imaginary 
time, the tunnelling barrier is created by the outgoing particle itself, whose 
trajectory is from the inside of the black hole to the outside, a classically for- 
bidden process. The probability of tunnelling is proportional to the exponential 
of minus twice the imaginary part of the action for this process in the WKB 
limit. Because of energy conservation, the radius of the black hole shrinks as a 
function of the energy of the outgoing particle; in this sense the particle creates 
its own tunnelling barrier. 

In the spherically symmetric case the emitted particle is taken to be in an 
outgoing s-wave mode, and so we use the plus sign in ((HJ. At the horizon (where 
g(r) = f{r) = 0) then r = provided is well defined there. The imaginary 
part of the action for an outgoing s-wave from ri„ to ro^t is expressed as 

Im / = Im / prdr = Im / / dp'^dr (7) 
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where Tin and Tout are the respective initial and final radii of the black hole. 
The trajectory between these two radii is the barrier the particle must tunnel 
through. 

We assume that the emitted s-wave has energy uj' « M and that the total 
energy of the space-time was originally M . Invoking conservation of energy, 
to this approximation the s-wave moves in a background spacetime of energy 
M ^ M — Lo' . In order to evaluate the integral, we employ Hamilton's equation 
f = j^\r to switch the integration variable from momentum to energy {dpr = 
^), giving 

'"''^^ dr ,„ r r--' dr 



/ = / / -dH= / / -(-de.') (8) 
Jn„ Jm r Jo Jr^^ r 

where dH = —doj' because total energy H — M — oj' with M constant. Note 
that f is implicitly a function oi M — uj'. For the special cases where this 
function is known (eg. Schwarzschild) the integral in eq. ^ can be solved 
exactly in terms of a; |^. Another generalization of the null geodesic method 
spacetimes with a well defined ADM mass are considered (since dependence 
oi M — oj' is explicitly known) in order to obtain self gravitation effects; for our 
considerations self gravitation will be ignored^ . 

In general we can always perform a series expansion in uj in order to find the 
temperature. To first order this gives 

" dr , . „ dr ^, ^ 



(9) 



^0 Jr,„ r(r, M-cj') J^^^ r[r,M) 



OJ 



dr 



To proceed further we will need to estimate the last integral. First we note that 
Tin > rout because the black hole decreases in mass as the s-wave is emitted; 
consequently the radius of the event horizon decreases. We therefore write 
rin = ro{M) — e and rout = ro{M — uj) + e where ro(Af) denotes the location 
of the event horizon of the original background space-time before the emission 
of particles. Henceforth the notation tq will be used to denote ro(A/). Note 
that with this generalization no explicit knowledge of the total energy or mass 
is required since rg is simply the radius of the event horizon before any particles 
are emitted. 

There is a pole at the horizon where f = 0. For a non-extremal black hole 
/'(ro) and g' {r^) are both finite and non-zero at the horizon, so for these cases 
i only has a simple pole at the horizon with a residue of , ^ Hence 

\/ f {1-0)3' {ro) 

the imaginary part of the action will be 

Im/=^= +0{uj^) (10) 



^ See ref . 1151 for a discussion of self-gravitation effects in the context of the information-loss 
problem. 
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Therefore the tunneUing probabiHty is 



r = exp (-2Ini/) = exp (-/3w) (11) 
and so the Hawking temperature Th = /3~^is 



rr. Vf'{ro)9'{ro) , 
= (12) 

It is easy to confirm that for a Schwarzschild black hole the correct result 
of Th = gjj^^ follows. Situations in which the horizon does not have a simple 
pole correspond to extremal black holes, and need to be handled separately. 
One conceptual issue that arises when applying either the Hamilton- Jacobi or 
null geodesic methods to the extremal case is due to the fact that the model is 
dynamic, so emission of a neutral particle from the black hole implies a naked 
singularity, in violation of cosmic censorship. We will discuss the extremal case 
in Section 4. 

2.2 Hamilton-Jacobi Ansatz 

We next consider an alternate method for calculating the imaginary part of the 
action making use of the Hamilton-Jacobi equation We assume that the 

action of the outgoing particle is given by the classical action / that satisfies 
the relativistic Hamilton-Jacobi equation 

+ = (13) 

To leading order in the energy we can neglect the effects of the self-gravitation 
of the particle. 

For a metric of the form 

ds^ = -f{r)dt^ + -^ + C{r)h,jdx'dx^ (14) 

the Hamilton-Jacobi equation H13|) is 

f{k ^ ^("^(^''^)' + wf'^^'^ ^ " " ^^^^ 

There exists a solution of the form 

I = -Et + W{r) + J{x') (16) 

where 

dtI = -E, drl^W'ir), dj^,h 
and that the J^'s are constant. Solving for W{r) yields 
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(for an outgoing particle) and the imaginary part of the action can only come 
from the pole at the horizon. It is important to parameterize in terms of the 
proper spatial distance in order to get the correct result [Tl] . 

For the first method the Painleve coordinate r was the proper spatial dis- 
tance. In this case the proper spatial distance between any two points at some 
fixed t is given by 

dff^ ^ — + C(r)h,,dx'dx^ (18) 
dir) 

As with the null geodesic method we are only concerned with radial rays, and 
so the only proper spatial distance we are concerned with is radial 

da' ^ 



Employing the near horizon approximation 



/(r) = /'(ro)(r-ro) + ... (19) 
g(r) = g'(ro)(r - ro) + ... 



we find that 



dr .,2^^ (20) 



V5W Va'iro) 
is the proper radial distance. So for particles emitted radially 




W{a) = " / -JE^ - ^g'(ro)/'(ro) (m^ + 

2TTiE 



(21) 



and from this point the computation is the same as for the previous method, 
yielding 

Th = ^ (22) 

for the temperature. 



3 Applications 

3.1 Rindler Space 

We first illustrate how these methods apply for the horizon of an accelerated 
observer. We shall employ different coordinate systems for 2D Rindler space 
to show that the same temperature results from applying the two tunneling 
methods directly. 

The forms of the Rindler metric being used are: 
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ds^ = -{a'x'-l)dt^+ dx^ (23) 

a^x^ — 1 

ds^ = -a^x^dt^ + dx^ (24) 

where a is the proper acceleration of the hyperboUc observer. Here there is 
no well defined total mass or energy as with Schwarzschild, but there are well 
defined horizons. The metric H23|) locates the horizon at a; = i, whereas for 
the metric (j24|l it is at a; = . 

We consider a null particle to be emitted from the Rindler horizon, and it 
is reasonable to assume the emitted particle will have a Hamiltonian associated 
with it. However providing an explicit definition for the total energy of the 
space-time is less than clear, though it has been claimed recently 27 that one 
can associate a surface energy density ct = ^ with a Rindler horizon and a 
total energy E = with the spacetime. In the context of the null geodesic 
method we expect that the Hamiltonian of the space-time will correspond 
to the total energy E (perhaps with respect to some reference energy via a 
limiting procedure) so as long as the emitted particles have w << in which 
case the method is applicable. We shall proceed under the assumption that 
we can use Hamilton's equation and follow through the derivation for the null 
geodesic method as before. We shall find that these assumptions are justified 
a-posteriori. 

The null geodesies for H24I) in the x-direction are given by 



X — ±ax 

and so 

r-'* dx TTUJ 

Im I = LO — = 

J^.^ ax a 

yielding a temperature of 

a an 

where the surface gravity at the horizon is k = 1 28 . 

We now employ the Hamilton- Jacobi ansatz^ for the Rindler metric 123|l . 

Here / = a?x'^ — l^g — and at the horizon /'(I) = g'{l) = 2a so 

using (|2I} 

a 

again giving a temperature of Th — 

We see that we can recover the expected value for the temperature of Rindler 
space given our assumptions. This could perhaps be regarded further evidence 
that a total energy E = can be associated with Rindler space. 

^For earlier work in the Rindler context along these lines see ref. llljl . 
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3.2 Charged-Kerr Black Hole 

Wc consider next the Kerr-Newman solution. The Kerr-Newman metric and 
vector potential are given by 

ds^ = -fir, 6)dt^ + — — - 2i?(r, 6)dtd(t> + Kir, e)d^^ + S(r, 9)de^ 

9{r, 0) 

Aa = -^[{dt)a - a" sin^ e{d(^)a] (25) 



f{r,e) 

9{r,0) 
H{r,6) 



A(r) - a2 sin^ 9 
A(r) 



S(r,0)' 

asm^eir'^ + - A(r)) 



E(r, e)=r^ + 0? cos^ Q 
A(r) = + + - 2Mr 

We assume a non-extremal black hole so that M ^ > + so that there are 
two horizons at r± = M ± \J — — e^. 

There is a technical issue in applying these methods because the metric 
functions depend on the angle Q. In order to account for this we can no longer 
just look a generic spherical wave; instead we will examine rings of emitted 
photons for arbitrary fixed 9 = Oq. In the end we will discover our temperature 
is independent of (as it should be). 

A naive first attempt utilizing the null geodesic method would be to consider 
the transformation 

dt = dT-, ^J'}'\^"]^ dr 



g{r,0o)fir,eo) 



This gives the equation 



ds^ = 0o)dT^ + 2^/f{r,eo)J - MrdT + dr^ 



1 



- 2Hd(j){dT - V ^ ^""'^"^ dr) + Kd(j)^ (26) 
\/f{r, 9o) 



whose radial null geodesies correspond to 
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There remain divergences in the dtdr and drdcj) terms at the horizon, and 
is not weh behaved there. Only for sinfln = are these ehminated. 

Restricting further the calculation to 0o = or tt (in which case ^ = 1), the 
outgoing radial null geodesies along the z axis are 



r = l-Vl-g(r,0o)|sineo=o (28) 

which yields 

, dr 27rcj ^ + 

I = u) — = —-, rrTi = ZTTo;- 



f .g'(r+, 00 ) I sin e=o 2(r+-M) 
for the imaginary part of the action. This in turn results in the temperature 

^ 27rr2+a2 27r 2Af(Af + (AP a2 _ g2)i) _ g2 

which is the same as the found for the Kerr-Newman black hole by other means. 

The restriction to two specific values of 6*0 is because of the presence of the 
ergosphere. The calculation breaks down because /(r, ff) is actually negative 
elsewhere at the horizon (i.e. inside the ergosphere) and 9t is not properly 
timelike there . The two values = or tt correspond to where the event 
horizon and ergosphere coincide. 

To address this issue, we note that the original charged Kerr metric can be 
rewritten as 

dr^ H(r 0) 

ds^ = -F{r, e)de + — — + K{r, 9){d(j) - \'. dtf + ^ir)d0^ (30) 



F{r,e)^ f{r,0) + 



g{r,9) V ' Kir, 9) 

H^{r,9) A(r)S(r,( 



K{r,9) [r^ + a^)^ - A{r)a^ sinH 
where at the horizon 

H{r+,9) _ a 
K{r+,9) r\ + a? " 
So the metric near the horizon for fixed 6 — 6o 

ds^ = -FAr+,9o){r-r+)df+ ^''^ +K{r+,9,m-§f^,dtf 

9r{r+,9o){r -r+) K{r+,9o) 

(31) 



and defining dx^d(j)~ ^(r^M 

ds^ = ~Fr{r+,9o){r ~ r+)dt' + — — ^ + K{r+,9o){dx)^ (32) 

gr[r+,9o){r - r+) 

The metric (I32II is well-behaved for all 9o and is of the same form as Q with 
f{r) = Fr(r+,9o){r — r^) and ^(r) — gr{r^,9o){r — r^) . Hence we easily obtain 
the final result (fT^ 

^Frir+,9o)gr{r+,9o) 
47r 
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Explicit calculation of i^r(''+,^'o) and gr{r+,9o) yields 
A,(r+) 2r+-2M 

9r{r+,0o) - 



Frir+,eo) 



A,(r+)S(r+,go) ^ (2r+ ^ 2Af )(r^ + cos^(0o)) 
Although i^rl^^+j^o) and gr{r+,6o) each depend on Oq, their product 

is independent of this quantity. Hence the temperature is 
1 r+-M 1 (M2 



H 



2tt r\ + a2 27r 2Af (M + (Af2 - - e^) 5 ) - 



for any angle. 

We turn next to the Hamilton- Jacobi method to find the temperature. The 
action is assumed to be of the form 

I=-Et + J(t) + W{r, 0o) 

and rewriting this in terms of ^ 4' ^ flnt we find 

I=-{E- nHJ)t + Jx + Wir, 9o) 

where it is assumed that E — fin J > 0. This demonstrates a nuance overlooked 
in the null geodesic method; the transformation to x implies that E should 
be replaced by ii^ — ^hJ for the emitted particle. The reason for this is the 
presence of the ergosphere. The Killing field that is timelike everywhere is 
X = dt + ^nd^. A particle can escape to infinity only if PaX° < , and so 
—E + U,hJ < where E and J are the energy and angular momentum of the 
particle. 

Employing the metric in the near horizon form H31() . the final result for 
W{r, 6o) is the same as ()21|l with E replaced by E— fin J- 

again yielding the temperature over the full surface of the Black Hole 

1 r+ -M 1 [A-P -a^ -e^)^ 



T, 



H 



27r r\ + a2 2tt 2M{M + {A'P - - e^)^) - 



in full agreement with the previous method and with Euclidean space tech- 
niques. 
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3.3 Taub-NUT-AdS 

The general Taub-NUT-AdS solutions with cosmological constant A = —i/P 
are given by |^ 



6 dr^ 

ds^ = -F{r){dt + An\fl{-)d^f + __ + (^2 ^ n^){^de^ + fl[6)d^^) (34) 
where 

, r^-n^ _2Mr-fi(r4 + 6nV-3n4) 

^(^) = 2_L 2 + 2_L 2 (35) 

r -|- n r -\- n 

and fc is a discrete parameter that takes the values 1,0,-1 and defines the form 
of the function fk{d) 

sin 9 for k = 1 
fk{e)={ 6 fork = (36) 
sinh6' for k = —1 

One of the interesting properties of Taub-NUT spaces is the existence of 
closed timclike curves (CTCs) T7|. For these cases it is not clear how to apply 
the null geodesic method, since the emission of an s-wave particle would have 
to recur in a manner consistent with the presence of CTCs. 

However there exists a special subclass of Hyperbolic Taub-NUT solutions 
(for An^/f < 1 ) that do not contain CTCs. A discussion of Taub-NUT space 
and the special cases without CTCs appears in the appendix. We shall consider 
these cases in what follows. 

The temperature can be successfully calculated using the metric in the fol- 
lowing form: 

ds^ = -Hdt" + ^ + G{dip - ^^"/fc^i) rft)2 + (^2 ^ „2)^^2 (37) 

where: 

ic 2 j4/-eN 

H{r,e)^{F + F^'^:^^) (38) 
G(r, 6) = Af/-) [r^ + - fl{^-){An^F + k{r' + n^))^ (39) 

As before, we will consider rings at constant and use the near horizon ap- 
proximation. 



At the horizon 



/!(¥) [ ^4((r^+n2)cos2(f )), k^l 



4((r^ + n2)cosh2(%)), k = -l 
A{r\+n'^), fc = 
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Only when k — 1 (for which CTCs are present) and 6q — n (i.e. when cos(^) = 
0) are there any potential divergences at the horizon. Since 

Hrir+,0o) = F,{r+) 

the metric near the horizon for fixed 9 ~ Oq is 



ds^ = ~Fr{r+){r - r+)dt^ + 



Fr{r+){r - r+) 
G(r+,0o) 



G{r^M{d^ - ^-(;+)^"f'(^) (. - r^)dtf (40) 



= -F,(r+)(r - r+)de + -—^ ^ + G(r+, 0) V (41) 

Notice that defining x = ~ as with the charged Kerr case is pointless 
since VLh = 0. From this point the steps are the same as for the general 
procedures outlined for either the null-geodesic method or the Hamilton-Jacobi 
ansatz. Inserting this into the final result for temperature (either H12(l or 1)22(1 ) 
yields 



T„ ^ %1 (42) 

which is the same form found using the Wick rotation method |^ . 

To demonstrate this is straightforward. Consider the hyperbolic case (fc = 
— 1). The mass parameter can be written in terms of the other metric parameters 
upon recognition that F (r+) — yielding 

M = !^ '-^ - 

2Pr+ 

Using this mass in H42|l yields an expression for the hyperbolic Taub-NUT 
temperature of 

Comparing this to the result (|23]) for the hyperbolic Taub-NUT tempera- 
ture obtained from Wick rotation methods 

T,^^^ = ^ (44) 

(where N is the Wick rotated NUT charge) we obtain agreement upon recog- 
nizing that = — due to analytic continuation. Note however that there 
is an implicit analytic continuation in the definition of r+, since F{r^,n) 
F{r+,iN) (24|. 
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We close by commenting that although we considered only the fc = — 1 case 
to avoid problems with CTCs, both the A: = 0, 1 cases can be formally carried 
through, yielding the result (02) ■ In the context of the null geodesic method 
this situation could perhaps be interpreted by noting that Hawking radiation 
yields a thermal bath of particles, whose existence can statistically be reconciled 
with the presence of CTCs. In the context of the Hamilton-Jacobi ansatz the 
physical interpretation is less problematic provided the classical action for the 
particle can be considered to obey the Hamilton-Jacobi equation in the presence 
of CTCs. Our results suggest a-posteriori the answer is yes, but the matter 
merits further study. In this context we note recent work |29| demonstrating 
that there are no SU(2)-invariant (time-dependent) tensorial perturbations of 
asymptotically flat Lorentzian Taub-NUT space, calling into question the possi- 
bility that a physically sensible thermodynamics can be associated to Lorentzian 
Taub-NUT spaces without cosmological constant. Whether or not such results 
extend to Taub-NUT spaces without CTCs is an interesting question. 



4 Extremal Black Holes 

Extremal black holes need to be treated separately from the other generaliza- 
tions, since the integrand no longer has a single pole. The general results derived 
above are no longer valid and even the self gravitating terms may play a very 
important role. One of the properties that occurs in extremal case is the pres- 
ence of a divergent real component in the action. Although such a term does 
not contribute to the imaginary part of the action, this may be an indication 
that the tunnelling approach is breaking down and the calculation is becom- 
ing pathological. Unlike the Wick-rotation method, which involves finding an 
equilibrium temperature, the tunnelling approach describes a dynamical sys- 
tem. In this latter context when a black hole is extremal the possibility exists 
that an emitted neutral particle may cause the creation of a naked singularity, 
in violation of cosmic censorship. 

Such a pathological situation would be prevented if the tunnelling barrier 
had infinite height. However we do not find this to be the case, and an evaluation 
of the imaginary part of the action yields a finite temperature. This is consistent 
with the proposal that extremal black holes can be in thermal equilibrium at 
any temperature |3U| . 

For concreteness, we shall consider the particular case of the Reissner- Nordstrom 
metric , though we note that a diverging real component has also been seen to 
occur with the extremal CHS solution |14j. 

4.1 Extremal Reissner-Nordstrom 

The Reissner-Nordstrom space-time is described by the metric 
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The black hole is non-extremal when > and extremal when Q = M. 
For the non-extremal case when the tunnelling approach yields a temperature 
of 



. 1 yg^g^ 

using either of H12() or (|22|l . Note that the limit Q ^ M gives a temperature of 
zero. 

For the Reissner-Nordstrom case self gravitating effects have been calculated 
exactly [H] and the full emission rate is 



p ^ ^ ^-27r(2a;(M-f )-(Af-t^)V(j\/-t:j)2-Q2+M^M2-Q2j ^^^^ 

Expanding this emission rate in powers of to yields the temperature 1)46(1 to 
leading order. Note that setting Q — M yields a contradictory result, since 
the second term in the exponent becomes imaginary. This unphysical situation 
corresponds to an extremal black hole emitting a particle, a situation in violation 
of cosmic censorship. 

Consider a nearly extremal black hole that emits a particle so that the re- 
sulting black hole is extremal. This corresponds to substitution of Q = (M — u) 
where the black hole emits a null particle of energy u>. Insertion of this value of 
Q into (gnj yields 



1 VAf 2 - {M - 1 V2M^ 

Comparing this to the temperature obtained from the emission rate using 147|l 
gives 



r 



_ -27r(^2a)(A/-f ) + A/^A/2-(A/-(j)2) _ -27r(AfV2AnU+2wAf+0(wi )) 



From the definition Ijlll) we find that the temperature that is 0{y/uj) and again 
approaches zero the closer the original black hole is to extremality. Explicitly 



2nMV2M^ 4^ ^ ' 

which differs from the value given in 148(1 by a factor of 1/2. This discrepancy 
arises due to an inappropriate expansion implicitly used in obtaining ((48(1 , which 
assumes that uj « .^m ' invalid assumption for Q = [M — ui). In this 
context we note earlier work demonstrating that the transition probability of 
emitting such a particle that will make the black hole extremal is zero Q- 

We obtain a temperature that depends on the energy of the emitted particle. 
We pursue the extremal case further by considering a direct attempt to find the 
temperature from the metric in its extremal form 

ds' = -(1 - ^ fde + ^^^^ + r-dV? (50) 
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Using the Hamilton- Jacobi Ansatz as a first attempt at the calculation yields 
only a diverging real component, i.e. 



f{r) = 9{r) = Jpir- Mf + 0((r - Mf) 



so that 



dr 



= M 



dr 



[r - M) 



M ln(r - M) 



r — M = e M 



where M < r < 00 implies that bounds on a arc now — 00 < cr < cxd. Rather 
than considering an observer at infinity, we will consider an observer outside the 
horizon at some ri corresponding to <j{ri). From 117|) 



W{r) 



dr 



Af 



= M 



= M 



'"('■1) da 



1 



cr(ri) 



da, E^e 



1,2, h'^J.Jj, 



(51) 



For convenience we choose a{ri) so that the term under the root is never nega- 
tive. This integral is diverging and real, suggesting that no particles are emitted 
|14| . However this result is suspect in that it may be contingent on employing 
the near horizon approximation in the early stages of this method. 

We turn next to the null geodesic method. The outward radial null geodesic 
is given by 



r = l 



r 



1 



:{r-M)' 



1 

Insertion of this into @ yields 



Im / ~ Im 



{r - + Oiir - Mf) 



(52) 
(53) 



de'^de 



2M^e e {1 + jjee ) 



-2ujM^ Im 
Im 



2i 



'ee'» M(l-h^ee*«) 



)de 



0(i)+4Mc.[ln(e-'^ + ^)]|J 



AMlo Im 



In 



^ ^ M 
1 _ 2e 



(2n+ l)47rA/w 



(54) 
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where we have written (r — M) = — ee'^ and n is an integer. The first part of 
the integral is a real contribution of 0(i) that diverges as e ^ 0. It does not 
contribute to the imaginary part of the action. The imaginary part of the 
action leads to a non-zero finite temperature 

^ 8.Af(2n+l) 

for any integer n. The extremal temperature is quantized in units of the 
temperature of a Schwarzschild black hole! 

Note that this result depends crucially on the inclusion of the third order 
term, whose evaluation depends upon assumptions of the choice of Riemannian 
sheet. Had we expanded the integral for small e , we would have obtained a 
value for the temperature given by n = —1 in eq. (|^ . ie a negative temperature 
for the extremal black hole. 

Obtaining many (finite-valued) results for the temperature is reminiscent of 
the proposal that an extremal black hole can be in thermal equilibrium at any 
finite temperature |3(Jj . However we can see that these strange results arise due 
to an inappropriate use of the WKB approximation in the null geodesic method. 
Although writing (r — A/) = — ee** is consistent with the the assumptions = 
ro(M) — e and rout = ro{M — uj) + e (where ro(M) denotes the location of the 
event horizon of the original background space-time) for a non-extremal black 
hole, in fact the quantity rout does not exist, since the extremal black hole 
cannot retain an event horizon upon emitting any neutral quantum of energy - 
its only option for future evolution would appear to be that of evolving into a 
naked singularity, which cosmic censorship forbids. 

These results seem to imply that for black holes near extremality one must 
consider the full self-gravitating results, where the emitted particle drives the 
hole toward extremality. For an already extremal spacetime both methods yield 
a diverging real component in the action. This could be taken to imply that no 
particle can be emitted (since the alternative is creation of a naked singularity) . 

Based on the results of this calculation it would be interesting to consider 
the emission of a specific charged particle that would cause the black hole to 
go from one extremal black hole to another extremal black hole. In that case 
there would be well defined horizons before and after emission. 



5 Conclusions 

We have examined and compared the two different approaches to the tunnelling 
method for finding the black hole temperatures. Our results indicate that 
the method is particularly robust for non-extremal black holes, yielding results 
commensurate with other methods for Rindler space, rotating black holes, and 
Taub-Nut black holes. In this latter instance we have provided independent 
verification of the temperatures obtained for Taub-NUT spaces without CTCs 
via analytic continuation methods. Indeed it is not too difficult to show that 
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the temperatures even match when CTCs are present, though in this case an 
a-priori justification for the method is unclear. 

We also investigated extremal black holes, for which the tunnelling method 
is somewhat more problematic due to its dynamic nature. We found that the 
temperature is proportional to the energy of the emitted particles for black holes 
close to extremality. We also found that both methods yield a divergent real 
part to the action for extremal black holes, which is suggestive of a full sup- 
pression of particle emission. However the null geodesic method has a nonzero 
finite imaginary parts, whose values yields a countably infinite number of pos- 
sible finite temperatures for an extremal Reissner-Nordstrom black hole. This 
rather strange result arises because of a breakdown of the WKB method in the 
null geodesic approximation. This suggests limitations on the method, whose 
study would make an interesting subject for future work. An interesting test 
case would be that of emission of charged particles from an extremal black hole. 

6 Appendix 

6.1 CTCs and Taub-NUT space 

The presence of closed timelike curves in Taub-NUT space can be seen by 
considering the curve generated by the Killing vector and by examining g^p^ 

So for k = 1,0, and k = —1 with An^/d'^ > 1 the quantity g^^p < 0, yielding a 
timelike d^; the curve r = ro, t = to, and 6 = 6q becomes a CTC. 

However there is a range of hyperbolic Taub-NUT solutions that occur when 
Ati^/P < 1 that don't contain CTCs. Now it is possible for g^^ to be negative 
when An? /i"^ < Ibut this occurs for small values of ro and happens inside the 
horizon. Explicitly when k = —1 then g^^ is given by 

g^^ = 4sinh2(-)(r2 +n2) ( cosh2(-) - -^—-^ sinh2(-) j 

So g^^ > will always be true as long as < 1- Figures 1-3 are plots of 

1 — ^r^n'^ for a range of mass and NUT-charge . On the plots the x-axis is 
r/n. The k = —1 case corresponds to hyperbolic solutions whose event horizon 
has radius rt > n. Since g^^ only becomes negative when r < n (within 
At? I^? < 1) any CTCs are contained within the horizon (provided the mass is 
positive). So no CTCs are present outside of the horizon for the hyperbolic 
case when An'^/l'^ < 1. 
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n/L=0.25, M/n=0 

n/L=0.25, M/n=0.01 

n/L=0.25, M/n=0.1 

n/L=0.25, M/n=l 

Figure 1: Plots when 4n^/L^ < 1 for a range of masses 




n/L=0.5, M/n-0 

n/L=0.5, M/n=0.01 

n/L=O.S, M/n=0.1 

n/L-0.5, M/n=l 



Figure 2: Plots for when 4n^/L^ = 1 for range of masses. 
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n/L=0.05, M/n-0 . 1 

— n/L=0.5, M/n=0 . 1 
n/L=0.4, M/n=0 . 1 
n/L-0.25, M/n=0 . 1 

— n/L=0.10, M/n=0 . 1 

Figure 3: Plots for fixed mass and a range of n^/L^ 
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